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- . Abstract 

o 

Q \ Motivated by the existence of metal-insulator transition in one-dimensional 

non-interacting fermions in quasiperiodic and pseudorandom potentials, we 
studied interacting spinless fermion models using exact many-body Lanczos 

-. diagonalization techniques. Our main focus was to understand the effect of the 

O 

t^^ ' fermion-fermion interaction on the transport properties of aperiodic systems. 

We calculated the ground state energy and the Kohn charge stiffness Dc- Our 

Vh ' numerical results indicate that there exists a region in the interaction strength 

r^ ' parameter space where the system may behave differently from the metallic 

o' 

J-^ ' and insulating phases. This intermediate phase may be characterized by a 

Lj power law scaling of the charge stiffness constant in contrast to the localized 

C^ ' phase where D^ scales exponentially with the size of the system. 



I. INTRODUCTION 

Understanding the transport properties of interacting many-fermion disordered systems 
has been one of the most challenging problems in recent years. Non-perturbative effects 
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are introduced by strong fermion-fermion correlations occurring in these systems and the 
available numerical tools are plagued by the exponential increase of the Hilbert space as 
larger and larger systems are studied. 

One-dimensional (ID) interacting fermion systems have emerged as useful models since 
exact solutions exist for some of them. Additionally, non-interacting ID models, in aperiodic 
(such as quasiperiodic or random) potentials are known to exhibit interesting phase dia- 
grams. In case of random potentials, the non-interacting fermion systems exhibit Anderson 
localizationu. Anderson insulator is characterized by a gapless point spectrum. This phase 
should be distinguished from the one occurring in a Mott insulator which is characterized by 
a finite gap originated in fermion-fermion interaction effects. The quasiperiodic and deter- 
ministic aperiodic models exhibit metal-insulator transitionoQ. In addition, these aperiodic 
models exhibit a phase which is intermediate between the metallic and the insulating phase 
and has been christened as critical. Critical behavior is characterized by multifractal quan- 
tum states and energy spectrumErQ. The transport properties for this phase, for example, 
resistance, are known to be oscillatingpcl. Unlike the Anderson localized phase characterized 
by the exponential decay of the single particle fermion wave function, the critical states 
exhibit at most power law decay. 

The possible existence of a phase which is in-between metallic and insulating in an 
aperiodic many-body fermion systems is a completely open problem. In this paper we present 
some preliminary results regarding this interesting problem in spinless many-body fermion 
model. We compute the Kohn stiffness constant describing the ground state conductivity 
at zero temperature by exact diagonalization on finite chains of various sizes. The main 
question that we address here is how the metal-insulator transition and the critical phase of 
aperiodic non-interacting models are affected by the presence of many-fermion interaction. 

The significance of our studies can be also viewed from a different perspective. Recently, 
there have been many theoretical investigations of the effects of fermion-fermion interaction 
on the problem of persistent currents in mesoscopic ringsE^'EJ. These studies were moti- 
vated by the fact that the free fermion theory underestimates the magnitude of the observed 



persistent current compared to the observed experimental valueO. The studies on spin- 
less modelsll3llj showed that the repulsive interaction always decrease the amplitude of the 
current. As pointed out by Giamarchi et al0, both attractive and repulsive ground states 
have charge density fluctuations which are easily pinned by the disorder. However, in the 
case of attractive ground state, superconducting fluctuations screen the disorder resulting 
in the enhancement of current. Our numerical results on aperiodic interacting systems are 
consistent with the above picture. Novel aspect of our results is the existence of a peak in 
the stiffness constant at a characteristic value of the attractive interaction suggesting the 
existence of a new type of phase in strongly correlated disordered fermion systems. 

In section II, we describe the basic model under investigation and give a short description 
of the method utilized to compute the charge stiffness D^. Section III describes the results 
obtained from our numerical simulations. Finally, section IV contains our conclusions and 
a discussion of the possible implications of this research. 

II. MODEL SYSTEM AND CHARGE STIFFNESS CALCULATION 

We studied an interacting spinless fermion model on a ID ring in an aperiodic potential, 

N N N 

H = - Xl(4ci+1 + 4+lCi) + Vj2^i'^i+1 + J2^i^i- (2-1) 

The site dependent potential is chosen to be of the form, hi = Xcos{2nai'^). Here, A 
represents the strength of the potential and a is an irrational number which is chosen for 
convenience to be the Golden Mean ( 2"^ )- "^^^ parameter u, determines the nature of 
aperiodicity: for u = 1 the potential is quasiperiodic while for z^ > 3 it has been shown to 
generate pseudorandom terms. Therefore, this particular form of potential facilitates the 
study of both quasiperiodic as well as pseudorandom cases by varying u. 

For A = 0, the interacting spinless fermion problem could be mapped to the Heisenberg- 
Ising XXZ spin problemlij. This is an old problem that have been extensively studied and 
for which a closed Bethe's anzast solution existsEj. In the non-interacting limit {V = 0), the 



quasiperiodic case (z/ = 1 ) can be reduced to the famous Harper equationE3. The Harper 
equation exhibits a metal-insulator transition in one dimensiorod. At the onset of transition 
Ac = 1, the quantum states are neither extended nor localized but instead exhibit fractal 
characters and have been termed as critical. The spectrum contains an infinite number of 
gaps and is believed to be a Cantor set of zero measure. These interesting aspects of the 
wave function and the spectra have been shown to be reflected in the transport properties 
such as Launder resistencea'Q. 

For finite values of fermion interaction V , the problem is more complicated due to the 
many-body nature of the wave function. Unlike the non-interacting case, where the behav- 
ior of the system can be described by studying the single particle wave function and the 
associated eigenvalues, in many-body problem, one needs an alternative method to charac- 
terize the nature of many-body state . Recently, the Kohn stiffness constant D^ has been 
introduced to characterize the difference between the metallic and insulating phases as it 
gives a direct quantitative measure of the electronic conductivity of the systemli3. In this 
paper, we will use D^. to determine the nature of the phase of the aperiodic system. 

To compute the Kohn stiffness constant Dc , we assumed periodic boundary conditions 
for the fermion model described by equation ( |2.1| ). We are interested in the persistent 
current response to a vector potential of magnitude \A\ = ^ in the x direction, where $ 
is the flux threading the ID ring and N is the number of sites in the chain. We used the 
Lanczos diagonalization methodta to obtain the ground state energy of the system Eq{^) as 
a function of the flux. The Kohn stiffness constant is then given by the equationll30. 

The numerical calculations were done using several sizes for the ID ring up to maximum 
size of A^ = 14. We studied the behavior of the Kohn stiffness constant Dc as a function of 
the parameters V and A describing the strength of the interaction and the strength of the 
aperiodic potential respectively. 



III. SIMULATION RESULTS 

We did simulations of the aperiodic spinless system For various sizes A^, and the electronic 
densities p = jf-, for many values of the parameters V and A. To simulate golden mean 
quasiperiodicity into the model, we used Lanczos methods for systems of various Fibonacci 
sizes. Furthermore, in order to keep the fermion density p almost a constant, we worked 
with densities which are the rational approximants to the golden mean a or the square of 
the golden mean a^. This procedure provides several possible sizes (5, 8, 13) for which the 
Lanczos diagonalization can be done at almost constant density. Therefore, our studies are 
for systems away from half-filling where the umklapp processes become irrelevant and the 
system in absence of disorder is metallic. By studying few different sizes, we were able to 
monitor finite size effects that could be present. 

Figure 1, shows the results for D^ describing the interplay between aperiodicity (with 
1/ = 1) and fermion-fermion interaction. The ground state transport properties are obtained 
for both the repulsive as well as for the attractive many-body interaction. The values of 
Dc decrease as A increases for all values of the fermion interaction. Consistent with the 
previous results on disordered systems", the repulsive interaction is found to decreases the 
Dc while the attractive interaction is found to increase the stiffness. As seen from the figure, 
an interesting aspect of the attractive ground state is the existence of a very prominent peak 
around V = —2.5 that survives even in the regime where the non-interacting system is an 
insulator. The location of the peak is insensitive to the value of A, however the peak gets 
narrower with increase of A. Similar effects were also observed for other densities which are 
rational approximants of cr^. 

The figures indicate the possible existence of a region where Dc may takes intermediate 
values: between those of a metallic and those corresponding to the Anderson localized insu- 
lating phase. In order to study this effect in more detail, we performed further simulations 
to study the behavior of Dc versus A for different values of V (see figure Ic). The figure 
clearly shows that in spite of the general decreasing behavior of Dc with A for any value of 



V, for V = —2.5 the values of D^ are relatively larger than those for V = —1 and V = 0. 
This confirms the special behavior of Df. around V = —2.5. 

Figure 2 shows the corresponding results for the pseudorandom case {v > 3). For A = 
case we observe the well known transition from an insulating to the metallic phase as the 
interaction strength parameter V is varied. In analogy with the results of the quasiperiodic 
case, the attractive ground state in the pseudorandom case also exhibits the novel character- 
istics namely the existence of a peak in Dc as V is varied. However, unlike the quasiperiodic 
case, the effect is less dominant in the pseudorandom. For the repulsive ground state, the 
interaction decreases the Dc value. For larger values of A a slight increase is observed for 
small values of V. 

The intriguing behavior of the conductivity manifested by a characteristic peak as the 
attractive interaction V is varied, for various strengths of the deterministic disorder A was 
further analyzed by studying the dependence of Dc on the size of the system. We notice 
that even though the height of the peak in Dc decreases with A^, the size of the system, this 
decrease is rather slow, particularly in comparison with the variation in Dc with A^ in the 
regime far from the peak, i.e. in the insulating phase. We conjecture that in the insulating 
phase, the Dc decays exponentially with the size of the system, while in the regime near the 
peak, the charge stiffness decays as a power law. This conjecture was verified at a special 
point V = 0. For A = 2 where the system is known to be in a critical state, the Dc value 
was found to exhibit a power law decay with the size of the system. On the other hand, in 
the localized phase, the Dc value was shown to decay exponentially with A^. 

We would like to mention that the reasoning for the above conjecture regarding the 
variation in Dc with A^ is in our observation that figure 1 describing the dependence of 
Dc with V is reminiscent of the plots of total band width (TBW) in the models of non- 
interacting fermions in a quasiperiodic field. It is known, for example, that the TBW for a 
critical non-interacting fermion model decays like a power law with the size of the systems. 
This is in clear contrast to the behavior at the localized phase where the TBW decays 
exponentially. This fact offers a powerful criterion to distinguish both phases. Since the 
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TBW is a measure of transport properties, it is conceivable that its scaling properties are 
similar to those of the D^. constant. 

IV. CONCLUSIONS AND DISCUSSION 

The purpose of this work has been to gain some insight in understanding of the trans- 
port properties of a ID chain of spinless fermions under the concurrent presence of strong 
interaction and aperiodicity. The numerical Lanczos diagonalization method used for these 
simulations is in essence an exact methodtJ. Unfortunately, its major limitation is due 
to practical aspects; even for relatively small sizes of the systems, large amounts of RAM 
memory and fast CPUs are required to obtain accurate values of the ground state energy. 

In spite of these difficulties, we think that is possible to extract some valuable information 
from our numerical simulations as our results are independent of various parameters such as 
the amount of disorder, fermion density and the system sizes. First of all, our calculations 
confirm the known fact that the effect of disorder in the interacting system is the localiza- 
tion of the metallic phase. However, the simulations show that for aperiodic potentials, we 
may have a phase that is not completely localized as in an Anderson insulator, nor com- 
pletely metallic as in a Mott conductor. This fact seems to be independent of the density. 
Comparison of quasiperiodic and pseudorandom cases suggest that in the quasiperiodic case 
this intermediate phase may exist in a finite parameter interval while in the pseudorandom 
case this phase may exist at a single point marking the boundary between the metallic and 
Anderson insulating phase. This distinction between the behavior in the quasiperiodic and 
the pseudorandom cases is reminiscent of the analogous distinctions that is known to exist 
in some non-interacting models studied previouslyQ. 

Our preliminary results are on systems of rather small sizes, and therefore, it is hard 
to reach any definite conclusion regarding the the nature of the new phase proposed here. 
On the other hand, the existence of a region with a characteristic peak is seen on systems 
of various sizes as well as of various electron densities. The clarity with which this peak 



appears seems to hint some new mechanism involving some sort of competition other than 
the known screening of the disorder due to superconducting fluctuations. Therefore, it is 
tempting to speculate the possibility of a new type of phase in aperiodic strongly correlated 
systems. We hope that our preliminary studies will simulate further research in this area. 

ACKNOWLEDGMENTS 

The research of IIS is supported by a grant from National Science Foundation 
DMR 093296. JCC would like to thank the people at GIF (International Physics Center) 
in Bogota for their constant help and encouragement during his career and acknowledge the 
support obtained through a COLCIENCIAS-BID-ICETEX scholarship-loan from Colombia. 



REFERENCES 

^e-mail: jchaves@gmu.edu. 

^e-mail: isatija@sitar.gniu.edu. 

^E. Abrahams, P. W. Anderson, D. C. Licciardello and T. W. Ramakrishnan, Phys. Rev. 
Lett. 42, 673 (1979). 

4 J. B. Sokoloff, Phys. Rep. 126, 189 (1985). 

^H. Hiramoto and M. Kohmoto, Int. J. Mod. Phys. B 6, 281 (1992). 

6 1. I. Satija, Phys. Rev. B 48, 3511 (1993); Phys. Rev. B 49, 3391 (1994). 

^S. Das Sarma, S. He and X. C. Xie, Phys. Rev. B 41, 5544 (1990); J. C. Chaves, I. I. 
Satija and J. A. Ketoja, Physica Scripta 52, 614 (1995). 

8 Y. Liu and K. A. Chao, Phys. Rev. B 34, 5247 (1986). 

^ H. Cruz and S. Das Sarma, J. Phys. I France 3, 1515 (1993). 

^°G. Bouzerar, D. Poilblanc and G. Montambaux, Phys. Rev. B 49, 8258 (1994); M. Abra- 
ham and R. Berkovits, Phys. Rev. Lett. 70, 1509 (1993). 

"T. Giamarchi and B. S. Shastry, Phys. Rev. B 51, 10915 (1995). 

i^L.P. Levy, G. Dolan, J. Dunsmuir and H. Bouchiat, Phys. Rev. Lett. 64, 2074 (1990); V. 
Chandrasekhar, R.A. Webb, M.J. Brady, M.B. Ketchen W.J. Galager, and A. Kleinsasser, 
Phys. Rev. Lett. 67, 3578 (1991); D. Mailly, C. Chapelier, and A. Benoit, Phys. Rev. Lett. 
70, 2020 (1993); D. Mailly, C. Chapelier, and A. Benoit, Physica B, 197, 514 (1994). 

13 W. Apel, J. Phys. C 15, 1973 (1982); T. Giamarchi and H. J. Schulz, Phys. Rev. B 37, 
325 (1988); W. Lehr, Z. Phys. B 72, 65 (1988); R. Shankar, Int. Jour. Mod. Phys. B 4, 
2371 (1990); H. Pang, S. Liang and J. F. Annett, Phys. Rev. Lett. 71, 4377 (1993). 

-^^E. Fradkin, Field theories of condensed matter systems, Addison- Wesley Pub. Co., Red- 



wood City, Calif. 1991; A. Tsvelik, Quantum field theory in condensed matter physics, 
Cambridge University Press, New York, NY, USA 1995. 

^^ B. Sriram Shastry and B. Sutherland, Pliys. Rev. Lett. 65, 243 (1990); J. D. Johnson, J. 
Appl. Phys. 52, 1991 (1981). 

16 D. R. Hofstadter, Phys. Rev. B 14, 2239 (1976); P. G. Harper, Proc. Phys. Soc. London 
A 68, 874 (1955). 

i^D. J. Scalapino, S. R. White and S. Zhang, Phys. Rev. B 47, 7995 (1993); Q. P. Li and 
X. C. Xie, Phys. Rev. B 49, 8273 (1994). 

i^H. Q. Lin and J. E. Gubernatis, Comput. Phys. 7, 400 (1993); H. Q. Lin, Phys. Rev. B 
42, 6561 (1990); H. H. Roomany, H. W. Wyld and L. E. Holloway, Phys. Rev. D 21, 1557 
(1980). 



10 



FIGURES 
FIG. 1. Quasiperiodic case {u = 1). (a) Charge stiffness versus F for p = | and cr = |. (b) 
Cliarge stiffness versus V for p = -^ and a = -^. Different curves represent different values of A, 
namely: A = 0.0 (solid line), A = 0.5 (dotted line), A = 1.0 (short-dash line), A = 1.5 (long-dash 
line) and A = 2.0 (dot-dash line), (c) Charge stiffness versus A for p = j^ and a = -^. Different 
curves represent different values of V, namely: V = —2.5 (solid line), V = —1.0 (dotted line) and 
V = 0.0 (short-dash line). 

FIG. 2. Pseudorandom case {u > 3). (a) Charge stiffness versus V for /? = |. (b) Charge 
stiffness versus V for p = j^. Different curves represent different values of A, namely: A = 0.0 
(solid line), A = 0.5 (dotted line), A = 1.0 (short-dash line), A = 1.5 (long-dash line) and A = 2.0 
(dot-dash line). 
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